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Time-reversal symmetry breaking surface states in t — J model 
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Recently a phenomenological Ginzburg-Landau (GL) theory has been proposed to describe 
the occurrence of a locally time-reversal symmetry (T) breaking state near a Josephson junction 
between unconventional superconductors. In this paper we derive this type of GL free energy 
microscopically from the t — J model within a slave-boson mean-field approximation. The re- 
sulting GL free energy is shown to satisfy the conditions to have a T-violating surface state. 
The existence of this junction state may explain some of the recent experiments on High-T c 
superconductors. 
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§1. Introduction 

The symmetry of the superconducting state in high- 
temperature superconductors (HTSC) has been a sub- 
ject of intensive study as an important clue to clarify the 
mechanism of their superconductivity. The Josephson ef- 
fect allows us to investigate directly the phase properties 
of a superconducting order parameter (OP), and thus it 
is a powerful experimental probe for this study. Many 
experiments demonstrate that the superconducting OP 
in these systems has a predominantly d x 2_ y 2-wave char- 
acter, i.e. thfiJDP changes sign under 90°-rotation in the 
Cu0 2 planeBi). 

In ei-wave superconductors interface properties can be 
qualitatively different from those of conventional super- 
conductors because of the nontrivial angular dependence 
of their pair wave functions. We have shown that a 
locally time-reversal symmetry (T) breaking state can 
occur near an Josephson junctions between d-wave su- 
perconductofS and, in general, unconventional supercon- 
ductorstrErErtP. This T-violating state exists only near 
the surface and decays exponentially toward the bulk. It 
has important consequences on Josephson effects. The 
arguments which led to this conclusion were based on a 
phenomenological GinzbnmpLandaii (GL) theory includ- 
ing several assumptions. a H' Er E3' EiP 

In this paper we derive the GL free energy from the 
t — J model within a slave-boson mean-field approxi- 
mation, and demonstrate that it is possible to have a 
T-violating surface state. The reason we consider the 
t — J model is the following. Mean-field (MF) theories of 
the t — J model based on a slave-boson method predict 
a superconducting state with a d x 2^2^symmetryj!3 13 
and they ..may explain the magnetiaU'E£P as well as the 
transports properties of HTSC if the gauge fields rep- 
resenting the fluctuations around the MF solutions are 
properly taken into account. Thus, it is interesting to 
study whether the t — J model leads to a T-violating 
state, in particular, at the Josephson junction. 



§2. Mean Field theory and GL expansion of free 
energy 

We consider the t — J model on a square lattice with 
the Hamiltonian 



<i,j>a 



h.c. 



(2.1) 



<%,j> 



where the summation is taken over nearest-neighbor 
bonds and Cj CT = Cj CT (l — n^-o-)- We use the slave- 

boson method to enforce the condition of no double occu- 
pancy by introducing spinons (fi a ] fermion) and holons 
(pi; boson) (cjo- = b\fi a ). Then the Hamiltonian Js de- 
coupled with the following order parameters (0P)E3E3: 
(1) the bond OP, {b\bi) = xb and (f] a fi a ) = Xf which 
we assume to be homogeneous for all nearest-neighbor 
bonds: (2) the OP for the Bose condensation of holons, 
(bjbip; (3) the singlet RVB OP, - fjj^) = A£. 

The superconducting OP is given by the product of the 
last two, (bjbi)A*y In a slave-boson mean field theory 
there are four kinds of ordered states in all of which the 
bond OP are finite: (a) the uniform RVB state where 
only the bond OPs are finite; b) the spin gap state where 
the singlet RVB OP is also finiteS* ; In this state there is 
a (pseudo-) gap in the spin, but not in the charge excita- 
tions. Hence the name spin gap state; c) the supercon- 
ducting state where all three OP's listed above are finite; 
d) the Fermi liquid statar. A schematic phase diagram 
is shown in Fig.l. 

In this paper we consider only the optimally and over 
doped case where Tbe > Trvb, and so the critical tem- 
perature for superconductivity, T c , is given by Trvb- (In 
other words we do not treat the case where the onset of 
superconductivity is given by the Bose condensation of 
holons.) In this case we can take xb — ^ since we always 
consider the case T < T BE . Since the superexchangc 
interaction exists only for nearest-neighbor bonds, the 
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and the extended s-wave are natural candidates 



for the symmetry of the superconducting OP. The for- 
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mer (latter) is defined by A d (i) — (A^i +X — Ai^ +y )/2 
(A s (i) = (Ai,i +X + Ai : i +y )/2). Following a standard 
procedure we expand the free energy with respect to 
and A s .tZP The resulting GL energy is given after taking 
a continuum limit 



F 



d 3 



x 



[^{^■(^lA.f+^IA.^ + A-IDA,! 2 } 



]=d,s 



*2\ 



+ 7l |A d | 2 |A s | 2 + ± 72 (A3 2 A 2 + A 2 A 



+K{(D x A d )*(D x A s ) - {D y A d )*(D y A s ) 
+ C .c.} + ^(VxA) 2 ] 

(2.2) 

where D = V — i(2ir/$ )A with 4> (= hc/2e) being the 
standard flux quantum. The coefficients in F are given 

as 

3 J 

aj= T^-m- 



3 J x ^ tanhfo/2T) 

' E a w j <») 



7 = 
if, = 



(¥)sS'(&h 4 (*) 

^ 2 ^^^lsin 2 fc^ 2 (fc) 



(2.3) 



327V 
9J 2 
'327V 



E 



■ sin 2 k x u>d{k)ui s {k) 



where j = d or s, 71 = 2 7 , 72 = 7 /2, u d (k) = cosk x 
cos k y and ui s (k) = cosk x + cosk y . Here 



1 



1 



Wf 



X 2T' 



(2.4) 



and /(^fc) is the Fermi distribution function. 

The surface energy at the junction is calculated under 
the assumption of a specularly reflecting surface. We 
consider a planar interface parallel to the c-axis as shown 
in Fig. 2. In Fig. 2 the left and the right hand side are 
the same d-wetve superconductors described by the t — J 
model. Here the crystalline a-axis of the left hand side 
(L) is normal to the interface, while that in the right 
hand side (R) is taken as a free parameter, denoted as 
f (0 < f < 7r). (We can treat an S'/D-junction, where 
the left side is an s-wave superconductor, in a similar 
way. We consider this case in § 5.) In this configuration 
the important effects are associated mainly with the OP 
on the right hand side, so that we will simply represent 
the left hand side by a single ci-wave order parameter Ao 
only. 

The transmission and the reflection of electrons at the 
interface (/) may be described by the following Hamil- 
tonian, 



u ft,; 



(2.5) 



where ff£) (fu^ ) is the spinon operator for the left 
(right) side, and the matrix elements for tunneling (tk p ) 
and the reflection (r^p) are taken to be real. Treating Hi 
in a second-order perturbation theory we get the surface 
free energy Fj to lowest order in A's, 



Fj= Jds[ nM^s 

i,j = {d,s} 

+ Yl ^)(ASA i + A A*)]. 

i— {d,s} 



(2.6) 



The first term originates from the reflection of the 
Cooper pairs at the interface and the second term rep- 
resents the coupling between the two sides (gij = gji). 
For the £)/D-junction composed of the same supercon- 
ductors, the coefficients in eq.(2.6) are given as 

td= (if) J2 t kp Jl ^ k '^ UJd ^ LJd ^ 

kp 



/3J^~ 

VT 



kp 



9d 



V 4 



kp 



kp 



(2.7) 



+ (rt P +t 2 kp )J2^k,^ s (k) 2 



= (t 2 



.'/./.-= \ — ) E [dpMZk^p^sityuJdip) 

kp 



{rl P + 4p)M^k,QMk)uJ s (k)] 



with 



MM) 



, tanh( 



2T> 



tanh( 



2T> 



26 



2£ P 



tanh(|i) tanh(J§0 



(6+C P )(&"6) 2V & iv 

(2.8) 

The J\ terms represent the usual tunneling and the re- 
flection processes of a Cooper pair. On the other hand, 
the J2 terms give different types of reflection processes 
where one of the particle consisting of a Cooper is re- 
flected at the interface, while the other one tunnels to 
the opposite side (see in the Appendix for the derivation 
of the surface terms and the interpretation of J 2 ). 

Here we follow the method of Ref.l3E3 ) in taking the 
angular dependences of t kp and r kp , which are consistent 
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with the assumption of the specularly reflecting surface, 

v±(p) 



t 2 - 



P6(k n - p\\)6(k± - P± y 



-9(v ± (p)) 



v(p) 

r\ v = r 2 5(fc|| -p\\)5(k x + P± )^M9(v ± (p)) 

(2.9) 

where 9 is the step function. It should be noted here 
that if the angular dependences of the tunneling ma- 
trix elements were omitted, the coupling between two 
superconductors should vanish, and thus it would lead 
to unphysical results. 

§3. Relation to Phenomenological Theory 

Now let us briefly summarize the content of the phe- 
nomenological theory of Ref.4 to describe a T-violating 
surface state. The properties of the bulk system are de- 
scribed by the free energy F. We assume that only a 
d-wave OP is present in the bulk system, and the cou- 
pling between A d and A s is repulsive, if the latter were 
present. Namely, a d is negative below a critical temper- 
ature, while o s > for all T, and 71 —72 > and 72 > 0. 
An important point here is that positive 72 favors the 
relative phase between Ad and A s , denoted as tp ds , to 
be ±7r/2, i.e., the system would break T if both d and 
s- components coexist. We expect 72 > for the follow- 
ing reason. The complex combination d ± is, resulting 
from 72 > 0, would open a complete gap in the exci- 
tation spectrum, such that the system would gain more 
condensation energy. If (j) ds — or it {{d ± s)-state), the 
nodes remain though their locations arc shifted, and the 
gain of the condensation energy would be smaller. 

The interface properties are described by Fj. In the 
presence of an interface, A s can be induced by the reflec- 
tion of Cooper pairs and by the proximity effect from Ao 
in (L). These processes are possible because of the low- 
ering of the symmetry (translational and point-group) in 
the presence of the interface. In usual cases <p ds is deter- 
mined by the bilinear coupling terms in Fj which induces 
the finite A s . On the other hand, the 72 term favoring 
tfids = ±7r/2 is biquadratic and it would not be domi- 
nant in the sense of GL theory. Then cj) ds takes a value 
or 7r depending on the sign of td, t 8 and gd s - There 
is no T-breaking in this case. Under certain conditions, 
however, a different relative phase is favored leading to 
a state which breaks T. In Ref.4 it was argued based 
on the symmetry consideration that t d and gds should 
vanish at 9 = ir/4 and are small for 9 close to 7r/4, while 
t s remains always finite. Then, for 9 ~ tt/A <f> ds is deter- 
mined by 72, and, thus, it is possible to have a locally 
T-violating state near the interface. 

We numerically calculate coefficients in F and Fi as 
functions of the doping rate S and the temperature T. 
We take t/J = 3 throughout in this paper. First we 
solve the self-consistency equations for the uniform RVB 
state (the state where only xf and xb = 5 are finite), to 
get \F and the chemical potential, Xp- By using them 
we calculate the coefficients in F and Fj . For the doping 
rate where the superconductivity is observed in real sys- 
tems, i.e., < S < 0.3, only ad can be negative and all 
other coefficients of F are positive definiteEH>. This im- 



plies that we deal with a single d-wave component. We 
find indeed a positive 72 which favors the phase differ- 
ence <j>ds — ±7r/2. Hence, the above results are consistent 
with the assumptions of the previous purely phenomeno- 
logical theory. 

Next we consider the interface terms Fj. By taking 
the angular dependences of tk p and rk P as in eq.(2.9), we 
get the the coefficients in Fj as functions of tp, i.e., the 
angle between the interface and the crystal a-axis of the 
right hand side. We can explicitly examine the following 
properties using the expressions of t$ and <7y (i, j — d,s): 

t s (tp±n/2) = t s (tp), t d (tp±Tr/2) = -t d (ip) 

t s (p±ir) = t s (cp), t d ((f ± 7r) = t d (ip) (3.1) 
t s (-p) = t s (cp), td(-tp) = t d (<p) 

and 

g u (<p ± tt/2) = g u (</?), g ds (p ± tt/2) = -g ds {v) 
gu(<p ± 7r) = gu (</?), g ds (f ± 7r) = gds(<p) 
9ii{- i P)= 9n{v), 9ds(-¥>) = 9ds{v) 

(3.2) 

where i = d, s. We have calculated t% and gij numeri- 
cally, as shown in Fig. 3, with a choice of parameters 
8 = 0.15 and T = 0.8T C (T c - 0.108J). The important 
point here is that t s does not vanish for any value of tp. 
This property is robust for any doping S and tempera- 
ture T. Another point is that \t s \ is much smaller than 
\td\ except very near ip = n/4. This can be explained as 
follows. The integrand in the expression of t s has a factor 
(cosk x + cos k y ), while that of t d has (cosk x — cosk y ) 2 . 
Since we are treating the square lattice system, the fac- 
tor (cosk x + cosfc y ) is small everywhere on the Fermi 
surface, for the doping rate as small as S = 0.15, while 
the factor (cos/c^ — cosfcj,) 2 can be large there. Hence 
\t d \ can be much larger than \t s \. (For tp — 7r/4, t d has 
to vanish by symmetry. Then \t s \ can be larger than \t d \ 
if p is close enough to 7r/4.) 

§4. T-breaking state and Surface Current 

In this section we analyze the appearance of the T- 
breaking state. For this purpose we minimize the total 
free energy, F to t — F + Fj, with respect to the OP, A s 
and Ad, and the vector potential A. We use a coordinate 
system (x,y), with x (y) being perpendicular (parallel) 
to the interface. Under the transformation from (x, y) to 
(x,y), only K terms are changed, and it transforms to 

K[cos2p{(D x A d ){D x A s y - {D s A d )(D § A s )*} 

+ sm2p{(D x A d ){DyA s y + (DyA d ){D x A s )* }] . 

(4-1) 

Since there is no spatial variation of the OP along in- 
direction, the GL equations are formally written as 



dF „ dF 

7TTT = O x - 



dAy x ) ^d{d x Ay x )) 



(i = d,s) (4.2) 



and 



J X 



J 11 



dF 
dA x 

dF 



0, 



(4.3) 



dA y 4;t 
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where B = VsAg — VyA x and A x — A x cosip — 
A y sin tp, Ay = A x sin tp + A y cos tp. The boundary condi- 
tions at the interface are 

t OF dFi 



^(V £ A 4 (i)) dA^i) 



= (i = d,s) (4.4) 



and 



B(x = 0) = 0. 



(4.5) 



Now we analyze the instability to a T-violating state 
at the temperature T*. We consider a junctions as de- 
scribed in Fig. 2, and assume that on both sides the su- 
perconductors have the same properties, in particular, 
the same T ce j. At T = T c , Ad and Ao (d-wave) get finite, 
and A s is induced simultaneously near the interface due 
to the t s and gds term. For tp = it/ 4, td and gd s vanish, 
so only the 72term determines the relative phase <pd s - 
Since 72 (> 0) favors (f>d s = ±tt/2, T-breaking occurs. 
Thus for tp = 7r/4, we have T* — T c , namely, T-violation 
should occur at T c , the bulk superconducting transition 
temperature, irrespective of any other details. 

Numerical calculation for the GL equations confirms 
this argument. We have solved GL equations (4.2) and 
(4.3) under the boundary conditions, eq. (4.4) and (4.5). 
In Fig. 4 we show the spatial dependences of A^, A s and 
4>ds for 8 — 0.15, i.e., so-called optimum doping. The 
results for other values of S is qualitatively the same. 
For tp = it/4, 4>ds always takes a value tt/2 or — tt/2, 
once T becomes lower than T c . As tp moves away from 
tt/4, td and gd s become finite and compete with 72, which 
favors T-breaking states. Since t& grows very rapidly as 
tp moves away from tt/4, T* drastically decreases as a 
function of \tp — ir/4\ (see Fig. 5.). 

However, we should note here that at very low tem- 
perature the region of T-violation can be much larger 
because of the following reason. If T is lower than T cs , 
there is a chance to have finite A s without using tunnel- 
ing terms. In this case bothi^ and t s can be smaller than 
72, with keeping A s finiteElP. Unfortunately T cs derived 
from the t — J model is quite low (of the order of 10 -3 J 
corresponding to T cs /T c d ~ 10 -2 ), so the region T < T cs 
is not accessable by the GL theory. 

It has been shown that in a T-breaking slain ±he sur- 
face current along the interface can flowHBe The 
condition for the minimum free energy requires the van- 
ishing of the current normal to the junction (i.e., the first 
equation of (4.3)). This condition in turn leads to a fi- 
nite current along the junction which induces a magnetic 
field. In Fig. 6 the spatial distributions of the surface cur- 
rent and the magnetic field are shown. The extension of 
the current and the magnetic field are of the order of the 
penetration depth A. Note that no pet current is present, 
i.e. the integrated current vanishes 2 -!!), 

J diJ y (x) = ^- [B(x = 0) - B(x = L)]=0 (4.6) 

where L is the length of the system (L ^> A), and we 
have used eq.(4.5) . 



§5. Case of S/D-j unctions 

So far we have considered only D/D-junctions, where 
both sides are the same d-wave superconductors. The 
S'/D-junction (where the left side is replaced by a dif- 
ferent superconductor with isotropic s-wave symmetry) 
can be treated similarly. 

We consider a conventional superconductor for the left 
side of the junction. We take the band width, the chem- 
ical potential and the magnitude of the order parameter 
in (L) to be the same as in (R) for simplicity. Our aim 
here is to compare the qualitative features of D/D- and 
<Sy_D-junctions, but not the quantitative comparison. 

The surface GL energy can be calculated in a similar 
way as in § 2. The modified expressions are simply given 
by the replacement 



LO d (k) 



1 



(5.1) 



in eq.(2.7) for td and t s . For the iS/D-junction \t s \ is 
larger than \td\ (for tp not so far apart from tt/4), in 
contrast to the case of D/£>-junctions. (Fig. 7) However, 
the value of \t s \ in this case is much smaller than that 
of \td\ for £>/D-junction. This is due to the same reason 
as we have \td\ 3> \t s \ for the D/D-junction: the factor 
(cosk x + cos k y ) is always small on the Fermi surface, 
and so \t s \ is reduced even if the OPs of both sides have 
s-wave symmetry. 

The phase diagram of surface states for S'/D-junction 
is shown as a solid line in Fig. 5. In Fig. 8 the in- 
dependence of the relative phases between A^ and A s 
(Ao), 4>ds {4>o) is a l so shown. We find that in the phase 
diagram the region of T-breaking states is larger than for 
the £)/D-junction. This is the case even when we assume 
a rather small value of lAol/IA^ 1 ^. In S'/D-junction 
|A S | can be sizable due to large \t s \, and then the 72- 
term (fourth order) which favors T-violating states can 
compete with the second order terms (td and gds) even 
when tp is not so close to 7r/4. In the case of D/ in- 
junction, however, td and gds become dominant once tp 
moves away from tt/4, leading to a very limited region 
of T-violating states. 

§6. Conclusions 

We have discussed the S/D- and D/D-interface states 
using a GL-formulation derived from the slave-boson 
mean-field approximation for the t — J model. The T- 
violating interface state was found in both cases for cer- 
tain orientations in agreement with previous phenomcno- 
logical studies. 

In addition we found that the T-violating state is more 
likely to occur on an S/D- than on a D/D-interface. This 
distinction is mainly due to the fact that the Josephson 
coupling to a d-wave superconductor is generally weaker, 
because of the angular phase structure of the pair wave 
function which leads to destructive interference in the 
tunneling. Furthermore, the shape of the Fermi surface 
plays an important role, in particular, for the extended 
S'-wave state. 

Within the GL theory the range of angles tp where 
the T-violating state can appear is rather small for the 
D/D-junction. However, it should be noted that for tern- 
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peratures below T cs this region in the phase diagram 
could extend to a wider range. Obviously, the GL for- 
mulation does not allow us to access this temperature 
region. In our approach the ratio T cs /T c d is of the or- 
der 10 -2 so that temperatures below T cs are definitely 
far away from the range of validity of a GL theory. The 
low-temperature region can only be treated by methods 
based on Bogolyubov-de Gennes equations or quasiclas- 
sical theory which are considerably more complicated if 
one intends to include more realistic microscopic details. 

Previously a series of phenomena were discussed, in 
connection with broken time reversal symmetry on in- 
terfaces. A first example are flux lines on the interface 
which do enclose neither integer nor half-integer multi- 
ples of the stand^cLilux quanta, but some intermediate 
fractional fluxesatrEj' . Although the observation of frac- 
tional flux by Kirtley et alO is entirely compatible with 
our discussion here, it is not clear whether the data could 
not be explained in an alternative way. The main prob- 
lem for this kind of experiments lies in the requirement 
of comparatively long homogeneous interfaces, a condi- 
tion hard to satisfy with present technology. On the 
other hand, it was also discussed that phase slip effects 
on short interfaces could be used as a testEZP . So far no 
experimental data are available for this type of effect. Fi- 
nally an important aspect of the T- violating state is the 
presence of spontaneous currents. Their magnitude is 
small, however, and together with screening effects they 
would not lead to a net magnetization. Thus only a very 
sensitive probe with high spatial resolution, smaller or of 
order London penetration depth, would be sufficient to 
observe this effect. Until now the only method which has 
successfully observed the small magnetic fields induced 
by a T-violating state are prion spin rotation measure- 
ments in zero external fieldEf . 

It is important to notice, however, that also surfaces of 
d-wave superconductors can yield states with locally bro- 
ken time reversal symmetry. Recent experiments indi- 
cate that this type of state might be realized at low tem- 
peratures on YBCO surfaces with [llOj-orientatiorO. 
The evidence is given by the splitting of the zero-bias 
anomaly in the-./V r -characteristics as discussed by Fo- 
gelstrom et al.c^P. Clearly similar effects due to the re- 
arrangement of quasiparticle states are also expected in 
T-violating interfaces as discussed by Huck et al. and 
could possibly be tested by spectroscopy with a scanning 
tunneling microscopeEIP . 
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Appendix: Derivation of surface terms 

The surface free energy Fj in §2 is derived in the fol- 
lowing way. We describe the transmission and and the 
reflection of electrons at the interface by the Hamilto- 
nian Hi in eq.(5). Then the total Hamiltonian is given 



by H + Hj. Second order perturbation theory gives the 
excess energy due to Hj: 

AF = ~J^ dr(T T H I (r)H I (Q)) (Al) 

where (• • -)o denotes the average with respect to H . We 
substitute eq.(5) into this expression and decouple it in 
terms of Green's functions of spinons for both sides 

AF = AFi + AF 2 

AFi = ~E T E [tl P {G^(p,ie n )Gf 2 (k,ie n ) 

+ G^(p,ie n )G^(k,ie n )} +ii p G? 1 (p,ie n )G$ i (k, i 
AF 2 = -$>E [tlpiG^ipJe^G^i^) 

kp e T i 

+ G% 2 (p,ie n )G*(k,ie n )} 



+ ^{G^ip, ie n )GUk, i£n) + G&(p, ze„)G 2 fl 2 (fc, 

(A-2) 

where the Green's functions are defined by 



Gii (&! = — 

G 22 {k, itn) = 

Gf 2 (k, «£n) = — " 



el + e k + \A£\* 

el + e k + \A£\ 2 
At 



(A-3) 



ei+e k + \^£? 

G 2 i(k,ie n ) = Gf 2 (k,ie n )* 

with A = R or L, e n = 7rT(2n + 1), = 
(3J/4)(A d u; d (fc) + A s Lu s {k)) and A£ = (3J/4)A cj d (fc). 
Now we extract the lowest order (0(A 2 )) terms. In AFi 
there are terms of the form A L A R and A fl A fl in the 
numerator. The former results in the ^-terms, and the 
latter leads to the part of the g^-terms. We can also 
obtain the 0(A 2 ) terms from the denominators of Gn 
and G 22 in AF 2 . These terms are usually discarded in 
the discussion of the Josephson effect, since it does not 
depend on the phase difference of A's if both sides of the 
superconductors have only one component of the order 
parameter. In the present case, however, this term de- 
pends on the phase difference of A^ and A s , and thus 
it is equally important as the one from AF±. Neglecting 
terms independent of A, we perform the e„-summation 
in eq.(A2) to get the following expressions 

w= E J ^ k ,a P ) [ti P {{A L p yA* + a£(a«d 

kp 

+ r 2 p (A«)*A«] 
AF 2 = Yl MZUpMp +r 2 kp )\A*\ 2 

kp 

(A-4) 



6 



Kazuhiro Kuboki and Manfred Sigrist 



The expressions for U and gij (eq.(7) in §2) can be ob- 
tained from eq.(A4). For example, ga is given by the 
coefficients of |A d | 2 in AF = AFi + AF 2 . 

The meaning of the J2-terms is now obvious. Since we 
got these terms by picking up |A^| 2 in G(\ and (di- 
agonal components of the Green's functions), only one 
particle is transfered from (L) to (R) (or (R) to (L)) in 
this process. This means that one of the electrons con- 
sisting of a Cooper pair tunnels to the other side, while 
the other one is reflected when the Cooper is scattered at 
the interface. Hence the process including tk P (tunneling 
matrix element) can lead to the suppression {gd and g s ) 
and the interference (gds) of the superconducting order 
parameters. 
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Fig. 1. Schematic phase diagram of the t — J 
model within a slave-boson mean-field approximation. 
T x , Frvb and Tbe are the transition temperature for 
the bond order, singclt RVB order and the Bose conden- 
sation, respectively. 

Fig. 2. Interface between d— wave superconductors. 
The lines on both sides indicate the crystalline a-axis 
(the c-axis points out of the plane). 

Fig. 3. Numerical results of the coefficients in Fi 
for D/D-junction: (a) td and (b) g^ — d,s). Here 
S = 0.15, T = 0.8T C , P = 8.0 and f 2 = 0.2. 

Fig. 4 Spatial variation of the order parmeters and 
their relative phase. Here 5 = 0.15, P = 8.0 and f 2 = 0.2 
and ip — 7r/4. 

Fig. 5. Phase diagram of surface states in the plane 
of T and tp. Here 5 = 0.15, P = 8.0 and r 2 = 0.2. The 
line is for the D/D (S'/D)-junction. (The T-brcaking 
region for D/D-junction is invisible in this scale.) 



Fig. 



Distributions of the surface current J v 



(A/cm 2 ) and the local magnetic field B (G). Here S = 
0.15, T = 0.8T C , P = 8.0, r 2 = 0.2 and ip = tt/4. 

Fig. 7 Numerical results of the coefficients in Fj for 
SyD-junction. Here 5 = 0.15, T = 0.8T C , P = 8.0 and 
r 2 = 0.2. 

Fig. 8 The relative phases between and A s (A ), 
4*ds (M as a function of <p. 
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